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Abstract 

The Faddeev equations for the three-body bound state are solved directly as three-dimensional integral equa- 
tions without employing partial wave decomposition. Two-body forces of the Malfliet-Tjon type and simple spin 
independent genuine three-body forces are considered for the calculation of the three-body binding energy. 
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1. Introduction 

Three nucleon (3N) bound state calculations are 
traditionally carried out by solving Faddeev equa- 
tions on a partial wave basis. After truncation this 
leads to a set of a finite number of coupled equa- 
tions in two variables for the amplitude. The cal- 
culations are performed either in momentum space 
[1-3], in configuration space [4,5], or in a hybrid 
fashion using both spaces [6]. Though a few par- 
tial waves often provide qualitative insight, mod- 
ern three nucleon bound state calculations need 34 
or more different isospin, spin and orbital angular 
momentum combinations [3]. It appears therefore 
natural to avoid a partial wave representation com- 
pletely and work directly with vector variables [7] . 
This is especially the case, if one wants to consider 



Email address: elster@ohiou.edu (Ch. Elster ). 
1 Corresponding author 



genuine three- nucleon force (3NF) effects. The true 
advantage stems from the fact that a 3N calcula- 
tion is carried out in Jacobi variables, whereas a 
typical 3NF has the form of two consecutive meson- 
exchange propagators between e.g. particles 1 and 
2, then particles 2 and 3. Using vector variables, the 
necessary coordinate transformations are numeri- 
cally realized by interpolations, whereas in a par- 
tial wave based calculation a large amount of cou- 
pling coefficients has to be evaluated and stored [8] , 
requiring large storage and memory capabilities of 
a computer architecture. 



2. Bound State Equation with a 
Three-Body Force 

The Faddeev component describing the bound 
state of three identical particles interacting via 
pairwise forces as well as genuine three-body forces 
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Fig. 1. Two meson exchange three-nucleon force 

can be written as 

jji = GaUP^i + (1 + G ti)G W^(l + P)<4) h (1) 

where i = 1,2, 3. In the following, we choose i = 1 
without loss of generality. The operator ti stands 
for the two-body t-matrix in the subsystem (jk) 
summing up the pair interaction in this system. 
The quantity shown diagrammatically in 

Fig. 1, is defined via 

Vi 23 = W {1) + W {2) + W {3 \ (2) 

where W {i) is that part of the full 3N force Vu 3 
which is symmetric for the exchange of nucleons j 
and k (j ^ i ^ k). The decomposition of Eq. (2) 
is natural, e.g. for the tttt exchange 3N force which 
is present in all currently available 3N forces. The 
free 3N propagator is given by 



G - h — q , 

m 4m 



(3) 



where p and q are the standard Jacobi momenta 
p^K^-kfe) (4) 



let* 



where ijk = 123 and cyclic permutations thereof. 
The permutation operator P is given as P = 
P12P23 + -P13-P23, and the full 3N wave function is 
related to the Faddeev component by 

|*) = (1 + P)|V). (5) 

The solution of Eq. (1) with — is described 
in detail in Ref. [7] and shall not be discussed here. 
For our model calculations we also use Yukawa in- 
teractions of the Malfliet-Tjon type here, however, 
we modify the interaction with a cutoff function of 
a dipole type [9] . 

In order to develop the algorithm, we concen- 
trate on a model 3N force of scalar meson exchange 
character, which can be written as 



g 2 s /A-nm s {MeW) A s (MeV) a s 
5.0 305.8593 1000.0 -1.73 

Tabic 1 

The parameters of the three-nucleon force. 

F(Q 2 ) F(Q' 2 ) 



W^ 1 



1 a 



S 9 2 



(2tt) 6 m s ys Q 2 +m 2 s Q'^ + m 2 ' 



(G) 



where m s represents the mass of the exchanged me- 
son. The momenta Q and Q' are defined in Fig. 1 
and given by 



Q = p p' ^(q q') 
Q' = p p' + ^(q q') 



(7) 



The function F(Q 2 ) = [(A 2 - rn 2 )/(A 2 + Q 2 )] 
represents a cutoff function for large momenta Q 2 . 
The parameters of the 3NF used in the presented 
calculation are given in Table 1. For evaluating 
Eq. (1) we need to calculate matrix elements like 



<pq|GW (1) |*> 



(8) 



Considering the momentum dependence in the me- 
son propagators of W^', one can easily imagine 
that the efficiency of any algorithm will crucially 
depend on the choice of coordinate system(s) when 
carrying out the integrations over the intermedi- 
ate momenta. Our numerical evaluation of Eq. (8) 
is based on the realization that can be inter- 
preted as two independent interactions, first in the 
subsystem (12), then in (31). Explicitly, we write 

( Pq |i^«|*) = J d 3 P ' J A' i(pq|p'q'>2 

d 3 p' 



p((p'-p") 2 ) 



(p' — p") 2 + m 2 
J d 3 p"'J dy" 2 <p"qV"q'"> 3 

3 (p""q'». (9) 



3 ,„, F((p"'-p 



d^p' 



(p'" - p"") 2 + TO 2 



Here the subscripts 1,2,3 of the bra and ket vec- 
tors define the meaning of the related vectors p 
and q, namely the particle number i singled out 
by qi as indicated in Eq. (4). We would like to 
point out that in each integration over a piece of 
the 3NF we only integrate over three variables, 
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gl/ATT m Q (MeV) A a (McV) g^/Air m r (MeV) A r (MeV) 
-3.5775 330.2104 1500 9.4086 612.4801 1500 
Tabic 2 

The parameters of the two nucleon pair force. 

namely the magnitude of a momentum and two 
angle variables. It turns out that it is most favor- 
able to choose momenta q' (q'") parallel to the 
z-axis, so that no interpolation on the unknown 
function ^(p,q' ,p ■ q') needs to be performed. The 
calculation of the transformation from one sub- 
system to another, e.g. 2(p"q'|p"'q / ")3> requires 
three-dimensional interpolations. We employ cubic 
Hermite splines [8], which prove to be both accu- 
rate and fast. 



3. Computational Approach 

The discretized Faddeev equation for a bound 
state (neglecting spin degrees of freedom) is an in- 
tegral equation in three variables on a typical grid 
size of 65*65*42 (momentum magnitudes, p, q, and 
angle between the momentum vectors). A priori 
the multidimcnsionality of the integral equation 
to be solved requires more memory. However, on 
an MPP system the number of variables and thus 
the memory do not pose a computational problem, 
since a variable defining a specific dimension of the 
grid can be distributed over a number of proces- 
sors, leaving a lower dimensional grid on each pro- 
cessor. As such, our three-dimensional approach 
is ideally suited as MPP application, and we can 
achieve an almost perfect load balance in our runs. 

The eigenvalue equation for the bound state 
is solved iteratively by using Lanczo's type tech- 
niques, here the method of iterated orthogonal 
eigenvectors [2]. For a typical run ten orthogonal 
eigenvectors are calculated per energy. We need 
about 5 to 7 energy iterations to find the ground 
state energy. 

The calculation of the kernel of the integral equa- 
tion means evaluating the matrix elements given 
on the right-hand side of Eq. (1) on a fixed grid p, 
q, and angle x — cos(p • q). The two-body t-matrix 
(with the two-nucleon pair interaction as the driv- 
ing term) is obtained by solving a system of linear 
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Fig. 2. The magnitude (in fm 3 ) of the 3N bound state wave 
function ^(p,q,x)) for x = 1 calculated with the potential 
parameters listed in Tables 1 and 2. 

equations of the form A • x = b, where A is typi- 
cally a 2500*2500 matrix. This system is solved for 
about 60 different vectors b, distributed over corre- 
spondingly many processors. The integrations over 
the 3NF terms are also distributed over the same 
number of processors, which depends on the size 
of the q-grid. Details about the different grid sizes 
and the dependence of the numerical accuracy on 
them can be found in Rcf. [7] . 



4. Results and Discussion 

For our model calculation wc use Yukawa inter- 
actions of the Malfliet-Tjon [10] type 



V(p',p) 

E 



(p'-p)2+m2 \K 2 a + (p' - p) 



A 2 



(10) 



with a short-ranged repulsive and a long-ranged 
attractive piece. The parameters used in our cal- 
culation are given in Table 2. A calculation of the 
three-body binding energy with this pair force 
gives £ t =7.699 MeV. The parameters of the 3NF 
are then adjusted such that its inclusion in the cal- 
culation gives a binding energy E t —&.590 MeV, a 
value close to the measured one. The full 3N wave 
function is calculated from the solution of Eq. (1) 
using Eq. (5) and shown in Fig. 2 as a function 
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Fig. 3. Contour slices of the 3N wave function at 2.4 fm 3 . 
The solid line represents the wave function calculated with 
2N pair forces and 3N forces, the dashed line stands for 
the wave function calculated with 2N pair forces only. 
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Fig. 4. Same as Fig. 3, except that the contour slices of 
the 3N wave function are taken at 10~ 3 fm 3 . 

of the momenta p and q at the fixed angle x = 1. 
A coparison of the 3N wave functions calculated 
with the 2N pair forces alone and with the inclu- 
sion of the 3NF is shown in Figs. 3 and 4. Both 
figures depict slices through the wave functions at 
specific values. 

The 3NF in our model calculation is attractive, 
as expressed by the larger binding energy. One 
would expect that the system becomes slightly 
smaller and acquires more high momentum com- 
ponents. Since the wave functions are both nor- 
malized to one, an increase of high momentum 
components will be seen as a decrease in low mo- 
mentum components, as seen in Fig. 3. Once the 
wave function decreases by some orders of magni- 
tude, the differences disappear. 

In summary, an alternative approach to state- 
of-the-art three-nucleon bound state calculations, 
which are based on solving the Faddeev equations 
on a partial wave basis, is to work directly with mo- 
mentum vectors. We formulate the Faddeev equa- 



tions for identical particles as a function of vector 
Jacobi momenta, specifically the magnitudes of the 
momenta and the angle between them, for the case 
when not only pair forces act between the parti- 
cles but also genuine three-body forces. As model 
forces we concentrate on scalar forces, a superpo- 
sition of an attractive and repulsive Yukawa inter- 
action for the pair force and an attractive Yukawa 
interaction for the three-body force. We demon- 
strate the numerical feasibility and accuracy of the 
solution. We want to point out that the incorpora- 
tion of the 3N forces is less cumbersome in a three- 
dimensional approach, and the algorithm can be 
made quite efficient on parallel architectures. 
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